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, We discuss a few aspects of the phase diagram of two color lattice QCD: we investi- 

1 gate the long distance screening analyzing the behavior of the interquark potential 

f**^ ' at large distances; we present a first set of results for vector mesons and diquarks; 

' we note similarities and differences between features at high temperature and high 

baryon density. 

(N 
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1 Introduction 

In this note we report on a numerical investigation of the phase diagram of 
I two color QCD in the temperature-chemical potential plane. Iruthis theory, 

^ ■ the fermion determinant is positive at nonzero baryon density ui3 and lattice 

simulations are possible. nnn 

Several studies of this model have already appeared Waa and, in partic- 
k>( \ ular, the phase diagram in the temperature-pdeasity plane at nonzero bare 

j_j ■ quark mass has been studied in some detail Qu'EI. There are three different 

\ phases: a hadronic phase with confinement and spontaneous (approximate) 

chiral symmetry breaking, a superfluid phase with diquark condensation and 
a plasma phase. For a nonzero quark mass and a zero diquark source, which is 
our case, the superfluid line is a bonafide phase transition, the. other lines are 
either crossover or first order transitions. It has been shown Em that, for some 
values of /i, either {ijjip) or ('0V') are non-monotonic function of temperature. 
An earlier, sketchy proposal of the phase diagram Q suggests a non-monotonic 
behavior for the diquark condensate. Such behaviors are rather common in 
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condensed matter (e.g. superconductors) but unusual in particle physics: that 
adds to the interest of two color QCD. 

Our aim is to characterize the various phases and their associated phe- 
nomena by continuing the study of the spectrum of mesons and diquarks and 
by exploring in detail the behavior of purely gluonic observables. 

Here we will discuss only a few aspects of our study: the behavior of the 
Polyakov loop and its correlators, and a first set of results for the vector meson 
ancLdiquark spectrum. A complete presentation, which will include topologi- 
cal a and thermodynamics results, will appear in a future publication IllJ. 



2 The simulations 



Numerical simulations of two color QCD with eight continuunx. flavors have 
been carried out on a 14'^ x 6 lattice by an HMD algorithm 0, for several 
values of temperatures, masses and chemical potentials. We have considered 
chemical potentials between zero and the lattice saturation and temperatures 
ranging from T ~ to T > Tc. 

Typically, we have run for 0(10000) MD steps, with dt = .02 and taking 
measures every 40 steps. 

A few test runs were performed with dt = .04 and dt — .08 in order to 
check that discretization errors are small at dt — .02. The accuracy of the 
inversion of the Dirac operator has been tested by use of the lattice Ward 
identity. 



3 Color screening 



It is interesting to remark that for two color QCD the Polyakov loop P is a 
real quantity: hence, the Polyakov loop correlator defining the strength of the 
quark-quark interaction (PP(i?))/(P)^, is the same as the quark-antiquark 
one (PP^(i?))/|(P)p . This property remains true also at 7^ 0. 

Note however, that this is the only special feature of the Polyakov loop 
and its correlators in two color QCD: all of the discussion can be carried out 
for a generic number of color N^. For this reason, it might well be that some 
of the finding presented here carry over to real QCD. 

The confined phase in the quenched approximation is characterized 
by an interquark static potential which rises linearly at any distance 
limi^^oo V{R) oc aR, a being the string tension. 

When quark degrees of freedom start to be dynamical a new phenomenon 
takes place. The color flux tube between two static quarks can now break by 
pair production of dynamical quarks popping out of the vacuum and screening 
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Figure 1. Polyakov loop correlation for /3 = 1.5 and various chemical potentials: screening 
and string breaking grow rapidly around /ic- The horizontal lines are drawn in correspon- 
dence of (P)^, demonstrating the cluster property of the correlators. 



the two static sources (string breaking): QQ — > [Qq){qQ)- The potential does 
not rise indefinitely anymore but flattens out at sufficiently large distance Rq. 

String breaking-in the hadronic phase has been observed at a moderately 
high temperature til, and it has also been suggested that the reaLjparticles 
present in a dense system further favors the breaking of the string 113. 

In the following, we will study the behaviour of the interquark potential 
as a function of fx, and we will contrast its behaviour with that observed as a 
function of temperature. 

Consider the Polyakov loop at spatial coordinates n — (ni, 712, na) 

P{n) ^ {l/N,)Tr [] Uo{n,r) . (1) 

T = 

We average P{n) on spatial planes to build: 

P3in:i) = il/N/) P(.ni,n2,n3) (2) 

"1 ,»2 = 1 
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Figure 2. Polyakov loop, and spatial Polyakov loop, as a function of /3 (left) and fj, 



Pi and P2 are similarly defined. 

By use of Pi{ni) we finally arrive at the zero-th momentum projected 
correlation function of the Polyakov line 

, 3 

^^P^^^^^ = ^ E E iP^MPhn^ + m ■ (3) 
i—1 Ui — l 

In Fig. we show the results for the Polyakov loop correlators at various 
values of the chemical potential. At high chemical potential we observe the 
expected enhanced screening, since the correlation functions approach a con- 
stant value at large distance. The gap between the plateaux ai fi — 0.2 and 
/i = 0.4 in Fig. |l| suggests the passage to a deconfined phase. We have then a 
direct evidence of the effect of the chemical potential on the gauge fields. 

We note that the trend, with chemical potential closely resembles that 
observed with temperature Ej: in both cases we have signals of deconfinement. 

In Fig. 2 the dependence of the Polyakov loop on /i (left hand side) and 
on T (right hand side) is contrasted with that of the spatial Polyakov loop, 
superimposed in the same plots. Again we note the similarity between the 
behaviour at high temperature and that at high density ; in both cases, only 
the temporal loop is affected. 

These results provide a direct observation of long range screening in a 
dense gauge system, obtained froni-a first principle study of the model, fur- 
ther confirming other observationsQ. They are thus providing an important 
confirmation of the standard expectation of the structure of a dense medium. 
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Figure 3. Polyakov loop versus the chiral condensate 



A more detailed analysis can give us more information on the screening 
properties. We note that some theoretical guidance can also be offered by 
low energy effective string models which provide the following expression 
(simplified for the sake of the present discussion) for the point-to-point corre- 
lators: 

{PP\R)) = ae-""^^' + (7V,iV/)2/3e-2'-°-^°^* . (4) 

The above equation shows the expected exponential dependence on the 
breaking length i?o, hence on the mass scale, and we also note the dynam- 
ical factors a and (3. It is then reasonable to assume that they depend on 
temperature and density: we will address this point in our future publication. 

Finally, in Fig. ^ we study the correlation between and {i^'ip) by 

plotting one versus the other. We plot there two sets of data: firstly, those 
obtained at fixed (3 = 1.5, and variable /x. Secondly, those at = 0, and vari- 
able (3. We see that chiral condensate and Polyakov loop are well correlated, 
for both data sets , and that the Polyakov loop looks consistently larger for 
the first data set ,/? = 1.5, /i ^ 0, supporting_the idea that a finite density of 
real particles further favors string breaking E3. 
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Figure 4. Vector mesons and diquarks propagators G{t) as a function of Euclidean time t. 
Upper, left: vector and pseudovector meson propagators in the hadronic phases. Vector 
diquarks are degenerate. Upper, right: as in left part of the diagram, but in the plasma 
phase: chiral symmetry is restored. Lower, left: vector mesons at high density; the propa- 
gators are degenerate, as it happens at high T, but chiral symmetry is still broken. Note 
anomalous spectral behaviour. Lower, right: vector diquarks at high density; time reversal 
symmetry is broken, and we note possible hints of vector condensation 



4 The vector sector 

The vector meson (qjiq) and diquark (q^iq) and antidiquark (qjiq) propa- 
gators at are constructed from the quark propagator, as explained in 
early work □ for the scalar spectrum, modulo the inclusion of the appropriate 
7 matrix (or, rather, of its staggered fermion representation). 

From the Pauli-Gursey symmetry of the SU(2) Action, the exact degen- 
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eracy of the vector (V) meson propagator with the pseudovector (PV) diquark 
propagator follows, as well as that of the pseudovector meson propagator with 
the vector diquank propagator. Analogous degeneracies have been noted in 
the scalar sector 

Finite density spectroscopy analysis has been already introduced and dis- 
cussed lJ. The following symmetries should hold true in the ensemble: 

G(t) = G{T - 1) for mesons (5) 
G{t, fi) = G(T - t, -fi) for diquarks (6) 

In addition to that, remember that staggered mesons contain states with dif- 
ferent lattice parity (the 'oscillating' channel): for instance, the vector meson 
contains both the p and the B particle, while the pseudovector contains the 
Al and again the p: 

Gv{i) - a(e(-"''*) + e'"''^^-*)) + (-l)*6(e(-™«*) + e"'^^^-*)) (7) 
Gpv{t) = a(e(""^i*) -I- e™-*i(^-*)) + (-l)*6(e(-'"''*) + e^'-C^"*)) (8) 

Since mesons do not 'feel' the effects of the chemical potential the masses 
of the forward and backward propagating meson states are identical. 

For /i ^ we expect different forward and backward masses in the diquark 
channels, reflecting the different backward and forward propagations in the 
dense medium. The occurrence of condensation, in either channel, should be 
signaled by the appearance of a constant term at large distance. 

Let us look now at the diagrams, Fig.^: the upper left diagram shows the 
vector and pseudovector meson propagators at low temperature, zero chemical 
potential. We have verified that they are degenerate with the vector diquarks 
of opposite parity. By increasing the temperature we observe chiral symmetry 
restoration in the vector channel (Fig.^, upper right diagram): all of the 
four vector states (mesons and diquarks) are now degenerate. At nonzero 
density mesons and diquarks with opposite parity are no longer degenerate, 
as anticipated. Let us follow independently mesons and diquarks then. 

In Fig.^ , left lower diagram, we see the vector and pseudovector mesons 
for jjL > fic- they are degenerate, as it happens at high T: of course this does 
not signal chiral symmetry restoration , since chiral symmetry is still broken by 
the diquark condensate. Simply, it tells us that the chiral condensate, which 
is responsible for the difference between these two propagators, is now zero. 
This is the same which happened with scalars. We observe also an important 
difference : the oscillating componenent is increased, as if the B particle 
contribution were more important at high density. This might indicate a 
significant difference of the spectral functions at high temperature and density. 
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suggesting that the spectral concentration of the p decreases at high density 
in comparison to that of the B meson. 

One last remark concerns the behavior of vector diquarks: in addition to 
the expected asymmetry, we observe (Fig.^ , right lower diagram) a behaviour 
which is not inconsistent with a possible plateau. Obviously, also a very 
low mass can mimick this behaviour, and all the usual caveat apply to this 
observation : first of all we should study colder lattices, where condensation 
phenomena are more prominent. Secondly, we should check finite size scaling. 
If this trend were to persist for a range fii lattice volumes then it could be 
interpreted as evidence for the predicted u vector condensation. 

5 Summary 

We have observed enhanced screening and string breaking in the critical region 
for chiral symmetry. We have shown that chiral condensate and screening 
properties are correlated both at finite temperature and finite density. These 
observations suggest that the (pseudo)critical line for chiral condensate and 
confinement run close to each other in the phase diagram of two color QCD. 
Moreover we noticed that, at comparable values of chiral condensate, string 
breaking is slightly more pronounced at finite density, and low temperature, 
than at zero density and finite temperature. This supports the viewJdiat 
recombination with real quarks in a dense system favors deconfinement E3. 

We have observed exact degeneracies in the vector sector of the spec- 
trum. We have studied the pattern of chiral symmetry for vector mesons and 
diquarks at finite temperature and density. Our results show a peculiar be- 
haviour in the meson vector sector at high density, noticeably different from 
that at high temperature: spectral concentration of the B particle seems dom- 
inating over the p, while at high temperature the opposite holds true. Finally 
we have noted possible hints of vector condensation in the diquark sector. 

One final remarks concerns the possibility to use two color QCD to test 
methods for simulating QCD at nonzero density. Sipme experiments with the 
Glasgow rewcighting have already been performedlij, while test of the analytic 
continuation from imaginary chemical potential are in progress cJ. 
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